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Abstract 

For Schwarzschild space-time, distributional expressions of energy-momentum 
densities and of scalar concomitants of the curvature tensors are examined for a 
class of coordinate systems which includes those of the Schwarzschild and of Kerr- 
Schild types as special cases. The energy-momentum density T ^ (x) of the grav- 
itational source and the gravitational energy-momentum pseudo-tensor density 
have the expressions T ^ (x) = —Mc 2 5®5q5^\x) and = 0, respectively. In ex- 
pressions of the curvature squares for this class of coordinate systems, there are 
terms like 5^(x)/r 3 and [<5*- 3 ^(a;)] 2 , as well as other terms, which are singular 



at x = 0. It is pointed out that the well-known expression RP a ^ v ' ({})R pa fMu({}) 
= 48G 2 M 2 /c 4 r 6 is not correct, if we define 1/r 6 = f lim e ^o 1/(V 2 + e 2 ) 3 . 
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1 Introduction 



In general relativity, many investigations have been made with regard to exact solutions 
of the Einstein equation and the singularity structure of space-time, but a distribution 
theoretical treatment of these has not been developed sufficiently This is the case even for 
the well-known Schwarzschild solution, which is given by, in the Schwarzschild coordinates 

ds 2 = - (l - (dx ) 2 + (l - t) {drf + f 2 [{d9) 2 + sin 2 W) 2 ] . (1-1) 

Here, a is the Schwarzschild radius a = 2GM/c 2 with G,M and c being the Newton 
gravitational constant, mass of the source, and the light velocity in vacuum Minkowski 
space-time, respectively. This metric, which is usually called the exterior Schwarzschild 
metric, is a solution of the Einstein equation in the vacuum region outside a static, spher- 
ically symmetric, extended gravitating body, and it is meaningless inside the body. 

The metric can also describe the gravitational field produced by a point-like 



particle located at x = 0. When we say, on the basis of the expression of the curvature 
square R pa ^ u ({})R papv ({}) obtained from the metric ( |1-1| ), that x = is a singularity of 
the Schwarzschild space-time, the source is considered to be point-like and this metric is 
regarded as meaningful everywhere in space-time. Then, we have 

T"(x) = , for x ^ , 



t M "(£) = 0, P M =(-Mc 2 , 0,0,0), (1-2) 



where T and t represent the energy-momentum densities of the gravity source and of 
the gravitational field, respectively, and P M is the total energy-momentum of the system. 
Thus, it is unambiguously expected that the energy-momentum density T of the gravity 
source has the expression 

T;(x) = -Mc\\^(x). (1-3) 

Nevertheless, this has not been given explicitly in any textbook of general relativity, as far 
as the present authors know. Also, in the literature, without giving distribution theoreti- 
cal examinations, the expression R pcrtJiU ({})R paflu ({}) = 12a 2 /r 6 is given, and the scalars 
R^{{})R^{{}) and R{{}) are not given explicitly, where R p ^x({}), ^(0) and R (U) 
stand for the Riemann-Christoffel curvature tensor, Ricci tensor and scalar curvature, re- 
spectively. 
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Recently, distributional expressions of energy-momentum densities of gravitational 
sources for the Schwarzschild and Kerr-Newmann space-times have been given by Balasin 
and Nachbagauerjl], |], || £| for the Kerr-Schild coordinate system. 

The purpose of this paper is to examine, from the distribution theoretical point of 
view, the energy-momentum densities, and the scalar concomitants R papu ({}) x R pafiu ({}), 

({})R pu ({}) and R({}) of the Schwarzschild space-time by assuming that the metric 
( |1 ■ 1| ) is meaningful everywhere in space-time. The examination is given for the coordinates 
f, 8, <p) and for the class of the coordinates r, 9, <fi) obtained by the transformation 

x° = tlx + f(r) , f = K(r) (1-4) 

with Q being a constant and / and K being functions of the coordinate r. (Note the 
following: (a)We obtain Kerr-Schild coordinates when f(r) = aln|a — r\,K(r) = r and 
Q = 1. (b)When f'(r) = 0,K(r) = r + a/2 and Q — 1, we obtain Lanczos coordinates 
which satisfy the harmonic condition.) 



2 Preliminary 

We briefly summarize the basics of general relativity, as a preliminary to latter discus- 
sion. 

In general relativity, the space-time is assumed to be a four-dimensional differentiable 
manifold endowed with the Lorentzian metric ds 2 = g^dx^dx" (/i, v = 0,1,2,3). At 
each point p of space-time, the metric can be diagonalized as ds 2 = rj^dX^p^dX^p 
with (r]^) = f diag (—1,1,1,1), by choosing the coordinate system {X M ; fi = 0,1,2,3} 
appropriately. 

The curvature tensor is given by 

^, {} )^„{;J-a,{;J + { A ;}{^}-{ A ;}{/J <»> 

with { p } being the Christoffel symbol. The fundamental action integral I is 

l=- c j (L G + h M )d 4 x, (2-2) 

where L^- is the Lagrangian density of a gravitational source and is the gravitational 
Lagrangian density given by 

L G d =^G. (2-3) 
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Here k is the Einstein gravitational constant k = 8ttG/c 4 and G is defined by 



clef 



a i r p i r a i r p 



with g = det(g^ u ). There exists the relation 

V=gR({}) = G + d^, (2-5) 

with 



A I „ A J /J 



Also, we have defined the scalar curvature byQ 

*({}) = (2-7) 

with the Ricci tensor 

R,A{}) = R X M{})- (2-8) 



From the action I, the Einstein equation 



<?/({}) = V(O) - \<R({)) = *V , (2-9) 



follows, where T^" is defined by 

T 

with 

T^mJpL (2-11) 

being the energy-momentum density of the gravity source. The energy-momentum pseudo- 
tensor density^ t of the gravitational field is defined by 

* M = <V L G - g aT ^ (2-12) 

OQ(TT,V 



with # CTriI , = f dg aT /dx". 



Raising and lowering the indices fi,v, A, • • ■ is accomplished with the aid of (<? Miy ) *= f (fly«/) 1 ancl i,9nv)- 
One may say, "It is meaningless to deal with f „ , because this is not a physically sensible local 
object." However, the following should be pointed out: (a)For an isolated system, t for an asymptotically 
Minkowskian coordinate system gives the total gravitational energy-momentum when integrated over a 
space-like surface. (b)Our main concern is not i„ , and it plays only an auxiliary role in our discussion 
(see also the footnote || (a) on page ||). 
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3 Schwarzschild coordinates 



Let us first consider the problem in Schwarzschild coordinates (x°,r,9,<f>) by which 
the Schwarzschild metric is expressed as Eq. (1-1). By using the Cartesian coordinates 
(a; , x 1 , x 2 , x 3 ), which are related to (x°, f, 9, 0) through the relation 

x 1 = f cos sin # , x 2 = f sin (f> sin 9 , x 3 = fcos9, (3-1) 

the metric takes the form 



ds 2 = g^dsfdx" , (3-2) 



where g^ v is given by 



fi'oo = -(1 - h) , g 0a = , 



g a p = S^ + hil-h)- 1 — , a, (3 = 1,2,3 (3-3) 



with h == a/r. 
We know that 



~ 


h! 


h 








f 




£; a 




i: 




= o, 






= 6J>( 


' h" 

v 2 



+ — Y --d> (3-4) 



G = 0, t„" = 0, (3.5) 

where the hatted symbols G, £ and represent, respectively, G, t and in the 
coordinate system {x M ; /x = 0, 1, 2, 3}. Also, we have defined h! = f dh/df and /i" = f d 2 h/df 2 . 
Regularizing the function h — a/r as a/ y/f 2 + e 2 with e being a real number, we obtain 

~o ae 2 
kT (x; e) = 



2(f 2 + e 2 ) 3 / 2 ' 



r 



o 



kT (x; e) = , KT a (x; e) = , 

~/3 , 3ae 2 i a i" ae 2 /5 e 2 \ 

K,T a {x;e) - d a - e ^ 5/2 - ^ 2 ^ 2 - ^ 5/2 ^- + — J , (3-6) 

where T (x; e) stands for the regularized T . Equations (3-4) and (3-5) have been ob- 
tained by the use of Eqs. (2-4), (2-12) and (2-9) with Eq. (2-10). 
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From Eq. (3-6), the desired expression 



T"(x) =\imT, V (x:e) 



follows, where use has been made of the relations^ 



X 



lim 



<y 2 ^ <y 2 j ^2^3/2 

r 2 



4tt(5 (3) ( 



x 



lim 



-2 _|_ e 2\5/2 



(3-7) 



lim 

e^O I f 



3 e 2 



+ 



2 (f 2 + e 2 ) 5 / 2 f 2 (f 2 + e 2 ) 3 / 2 



2vrrt (3) (a;) • 



(3-8) 



The regularized scalar quantities R({}; e), R^({}; e)R tlv ({}; e) and i?^({}; e) J R pCT ^({}; e) 
are calculated as 



#({};e) 

^({}; 6)^,(0; e) 
^({};e)4 ffMi/ ({};e) 



3ae^ 



+ 



2ae 2 



2 _)_ e 2\5/2 1 ^2^2 _|_ e 2^3/2 



3ae 2 



52 _|_ e 2p/2 



+ 2 



ae 



4a 2 



1 

+ T7 



^•2^i^j2 j ^2^3/2 

12a 2 e 2 



f 2 _|_ £ 2 |^^2_|_ e 2^2 j34j ^52 _j_ £ 2^4 

9a 2 e 4 



(f 2 + e 2 ) 5 ' 

respectively. Thus, the scalar curvature has the well-defined limit 

R({}) d =limi?({};e) = Ana6^(x) . 

e— >0 



(3-9) 



(3-10) 



However, the quadratic scalars do not have well-defined limits, which can be symbolically 
written as 



R^({})RM}) = \imR^{{} ] e)R, v {{} ] e)^A^ 2 a 2 [5^\x 



i^({})iW({}) = lim^({};e)^({};e) 

e— >0 



3 f 3 



(3-H) 



Proof of the relation (3-8) and a comment related to it are given in Appendix A. 
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with 1/f 6 = f lim e ^o l/{r 2 + e 2 ) 3 . Hence, we have 



^(0)^(0) + , R^{{})R P ^{{}) ^ 



12a 2 



(3-12) 



In Eq. (3-11), the terms (l/f 3 )5 (3) (&) and [S {3) ( a;)] are both ill-defined, and the symbol ~ 
is to denote the "equality" of the left- and right-hand sides in a rough sense. Mathematical 
formalism capable of describing these singular quantities is needed. 

It is worth mentioning here that the spherically symmetric non-scalar quantities 
f 3 i?^({}; e)R^({}; e) and r 3 Rf> a e)R pa ^({}; e) have the well-defined limits 



limf 3 F"({};e)^({}; £ ) 
limf 3 ^({};e)E w ({};6) 

e— >0 
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7ia 2 S^(x) 



12a' 



11 



-ixa~ 



» ( J(3)( 



X 



with 



def,. 
= hm 



e^o (r 2 + e 2 ) 3 



(3-13) 
(3-14) 



4 Generalization 



We now consider the problem in the coordinates r, 9, <fi) related to the Schwarzschild 
coordinates (x°, r, 6, 4>) through the relation (1-4). In the coordinates (x°, r, 8, 0), the metric 
has the form 



ds 2 = -A(dx ) 2 -2Ddx°dr+(B + C)(dr) 2 
+Br 2 [(d6) 2 + sin 2 6 (d^) 2 ] , 



(4-1) 



where we have defined 



A d =l f Q 2 (l-^- 
K 

a \ -i 



P 2 



C 



def 



■ - jf) <*T - £ - (i - £) if f • 



D S f! 



with[| p = f r, K' == dK/dr and /' == df /dr. The coordinates t = x°/c and r are time and 
space coordinates, respectively, only if 



def 



K, 

def 



(4-2) 



(4-3) 



4 The manipulation to denote r with p is for the regularization procedure employed below. 
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In the Cartesian coordinate system /i = 0,1,2,3} with x 1 = rcos^sin^, x 2 



rsin0sin#, x 3 = rcos#, the metric takes the form 



ds 2 = g^dx^dx" (4-4) 

with given by 



9oo = -A, g 0a = -D—, g ap = B5^ + C^. (4-5) 



We can obtain expressions for the quantities t , T ^ , R({}), -R^({})-R^({}) and 
i2 w "'({}).Rp 0/1I/ ({}) in terms of A, B, C, D, which are enumerated in Appendix B. 
By using Eqs. (4-2), (B-l) and (B-2), we can show that 



C = , (4-6) 



if Case 1: 

K(r) = Ar , (4-7) 

or Case 2: 

K(r) = Ar + a , f(r) = Tr + 5 (4-8) 

with A, T and S being real constants. (For Case 1, f(r) is arbitrary.) 

In what follows, we restrict our consideration to the above two cases, Case 1 and Case 
2 with positive^ A. In order to treat the singularity at r = in a distribution theoretical 
way, we employ the following regularization scheme: 

R.l First, replace K' with A in the expression of the function C . 

R.2 Then, replace K' and K" appearing in the process of calculation with Ap' and Ap", 
respectively. 

R.3 Finally, replace the function p with y/r 2 + e 2 . 

Then, we obtain the following results by using Eqs. (4-2), (B-l)~(B-4): 
Case 1 



_- IS 

5 Note the following: (a) It is expected that a delta function expression for T ^ (x) is obtained when 
= . (b) The variable r cannot be "radial" if A < 0. 
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(1) The regularized energy-momentum density T ^ (x\ e) is given by 



kT q (x; e) 



kT (x;e) 



nT a (x;e) 



j^i2 ^j^>2 j ^2^3/2 



0, KT a (x;e) = 



51 



3ae z 



ae 



2( r 2 + e 2)5/2 r 2 ( r 2 + e 2)5/2 r 



5 e 2 ' 

+ 



This leads to 



T."(x) =\imT"(x: e) 



p V") 



-Mc 2 \n\sX^( 



X 



which reduces to 



T:(x) = -Mc\\^(x) 



when \Q\ = 1. There is the relation 



T, (x;e) 



where we have defined 



Thus, we see that 



T» (x;e) 



T„ (x; e/A) 
<9£ A dx v p 



def 

c (9(x) cte^ dx p 



T x (x;e). 



^ (z) = lim 



e^O L 



For 



however, we have 



def <9(£) dx^dx^- p 

c ~ d(x)dx^dx^ X [X) 



T (x) 



T a (x) 



T?(x) 



-Mc 2 \Q\5 {3 \x) 



T (x) 



Mc 2 \n\ fl( ,x a . m , , 



in which the ill-defined quantity^ —6^ 3 '(x) appears, and hence 

T, V (x) ? 



T, (x) 



(4-9) 

(4-10) 
(4-11) 

(4-12) 

(4-13) 

(4-14) 
(4-15) 



6 Note that x a [(l/r)S^(x)] ^0= (l/r)[x a 5^(x)} and that £a=i <5 (3) (z) = 5 (3) (a:). 

the relation (A- 16) in Appendix A. 



(4-16) 

(4-17) 

See also 
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From Eqs. (4-14) and (4-17), we see the following: The correct expression for the 

— V 

energy-momentum density T ^ (x) is obtainable by transforming first the regularized 
density T ^ (£; e) and then taking the limit e — ► 0. However, it cannot be obtained, if 
the order of making coordinate transformation and taking the limit is exchanged. 

In Ref. 0, the expression T M (x) = -M<?8*8£8®(x) is given for Kerr-Schild co- 
ordinates. One might think that Eq. (1-3) can be easily obtained by transforming 
this to Schwarzschild coordinates. However, this is not the case; the situation here 
is similar to the one discussed above. 



(2) For the regularized scalar quantities, we have 

3ae 2 



R({};e) = A" 



+ 



2ae 2 



2 _|_ e 2)5/2 1 r 2( r 2 _|_ e 2)3/2 



R^({};e)R pu ({};e) = A" 



3ae 



2 l 2 



.2 _|_ £2)5/2 



+ 2 



ae 



R({};Ae) 

1 2' 



ly 1 ^ ^^«2 j ^2^3/2 



Ae)R p „({};Ae) 



(O;<0 = a- 



4a 2 



^2 j 



+ 



( r 2 _|_ £ 2"\2 1 r 4 



12a 2 e 2 
(r 2 + e 2 ) 4 



+ 



9a^e 



2^4 



(r 2 + e 2 ) 5 
^({};Ae)^({};Ae). 



(4-18) 



Equation (4-18) implies that the regularized scalars in the coordinate system /i 
= 0, 1, 2, 3} for Case 1 essentially agree with the corresponding scalars in Schwarzschild 
coordinates. 7 This is not self-evident, because we have made the regularization in 
each coordinate system independently.^ 



Case 2 



(1) The regularized energy-momentum density T ^ {x\ e) is given by 



~ o 
kT (x; e) 



\n\ 



6ae 2 



+ 



ae 



2 _|_ £ 2~)5/2 1 r 2^ r 2 _|_ £ 2)3/2 



7 Note that the product Ae in the hatted scalars in Eq. (4-18) plays essentially the same role as e, 

because the limit e — > is taken in the final stages. 
8 See also the paragraph at the end of this section. 
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KT a (x;e) 



\n\ 



a V") "~J 



n r 



-2 _|_ £ 2)5/2 



x a x /3 ae 2 



5 e 2 



p 3ae 2 

" a 2(r 2 + e 2 ) 5 / 2 + r 2 (r 2 + e 2 ) 5 / 2 \2 ' r 2 ) 



(4-19) 



which yields a limit having the same form as Eq. (4-10), upon using Eqs. (3-7), (3-8) 
and the relationP] 

. (4-20) 



lim 



ty ^^«2 j ^2^5/2 



The quantity 



~ is 

T (x; e) defined in the same way as Eq. (4-13) is given by 



T (x; e 



T (x; e) 



-\AQ\ 

, K 

|AO| 



ae 



r 2 [(Ar + a) 2 + e 2 ] 3 / 2 



Ar + a\ f 3a 



8i 



ae 



2 [(Ar + a) 2 + e 2 ] 5 / 2 
i^* 3 f 1 3 



[(Ar + a) 2 + e 2 ] 3 / 2 r 2 



(Ar + a) 2 2 (Ar + a) 2 + e 5 



(4-21) 



from which the limit 



lim 



T, (x;e) 



is obtained. Also, 



T, (x) 



= ^ T„ (x) (4-22) 
defined in the same way as for Case 1 has the expression 



T \x) 



T (x) 



T„ far 



-Mc 2 \AQ\ 



Ar + a 



5 {3) {Ax + a) 



Tf(x) 



-Mc 2 



\AQ\ ( Ar + a\ 2 x 
~0~ 



-5 (3 \Ax + a) 



(4-23) 



with a = ax/r, and we have a relation having the same form as Eq. (4-17). This 

— V 

and Eq. (4-22) imply that the correct expression for T ^ [x) is not obtainable for 



9 Proof of the relation (4-20) and a comment related to it are given in Appendix A. 
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Case 2 by applying the coordinate transformation and the limiting procedure to the 
regularized density T (x;e). 

The regularized scalar quantities are given by 



r 2 + e 2 



9ae 2 



2ae 2 



R^({};e)R, u ({};e) 



A(A V / ^T^ + a) 2 

(r 2 + e 2 ) 2 
A^Av^T^ + a) 4 1 2 

3ae 2 



( r 2 _|_ £ 2)5/2 r 2(j2 _|_ £ 2)3/2 



3ae 



+2 



+ 



(r 2 + e 2 ) 5 / 2 



ae 



( r 2 _|_ £ 2^5/2 1 r 2^ r 2 _|_ e 2)3/2 



R^({};e)R pa ^({};e) 



12a 2 



A 2 (A v / ^ 2 T7 2 + a) 6 
4a 2 

A 2 (A v / ^ 2 T^ 2 + a) 5 
2e 2 



A + 



ae 



+ e 2)3/2 



X 



+ 



A + 



9e 2 



ae 



(r 2 + e 2 ) 3 / 2 



r2v 92^2 ( r 2 + e 2 )3 /2 



A 2 (AVr 2 + e 2 + a) 4 
from which the following is known: 



4e 4 81e 4 
+ 



r 4 (r 2 + e 2 ) (r 2 + e 2 f 



(4-24) 



i2({}) d = f limit>({};e) = 0, 

e— >0 



^({};e)^({};e) 
i?^({};e)i?^({};e) 



2a 2 



A 2 (Ar + a) 4 r 4 (r 2 + e 2 ) 
12a 2 4a 2 



7^0, 



7^ 



(Ar + a) 6 A 2 (Ar + a) 4 r 4 (r 2 + e 2 ) 
12a 2 



(Ar + a) 



(4-25) 



Also, we can show that 



limr u r({};e)i? p ({};e) = 



limr^^({} ;e )V({}^) 



, for uj > 1 . 



(4-26) 



As is known from Eqs. (3-9) and (4-24), the regularized scalars in the coordinate 
system {x^; /z = 0,1,2,3} for Case 2 are not simply related to the corresponding 
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scalars in Schwarzschild coordinates. However, this is not a contradiction, because 
the regularization has been made in each coordinate system independently. 



5 Summary and comments 

In the above, we have examined the Schwarzschild space-time from a distribution 
theoretical point of view. Expressions of the energy-momentum densities T ^ and and 
of the scalars #({}), i^({})i^({}) and R papu {{})R P aM}) have been obtained for a 
class of coordinate systems including those of the Schwarzschild and of Kerr-Schild types 
as special cases. 

The results can be summarized as follows: 

(1) For Schwarzschild coordinates, employing the regularization scheme prescribed by the 

replacement a/r with a/\/f 2 + e 2 in the expression (3-3) of the components of the 
metric, we have shown the following: 

(l.A) The energy-momentum density T (x) of the gravitational source has the 
expression (1-3). 

(l.B) The scalar R({}) has the definite delta function expression (3-10), while 
R^({}) and RT» V ({})R pa M}) have the expression (3-11), which in- 
clude terms like [5^(x)] 2 and 5^(x)/r 3 . Also, r 3 i?^({}; e)^({}; e) and 
r 3 R pa ^({}; e)R pa ^({}; e) have the definite limits (3-13). The second relation in 
Eq. (3-12) shows that the well-known expression R pCT ^ v ({})Rp„^({}) = 12a 2 /f 6 
is not correct, as long as we define 1/f 6 = f lim e ^ l/(r 2 + e 2 ) 3 . 

(2) For the coordinates (x°,r, 0,0) related to the Schwarzschild coordinates (x°,f,#,0) 

through (1.4), we have shown the following: 

(2. A) If Case 1: K(r) = Ar, or Case 2: K(r) = Ar + a ,f(r) = Tr + E with 
A, T and S being real constants, the gravitational energy-momentum density 
vanishes, t„ = 0. 

(2.B) Restricting ourselves to cases with positive A and employing the regulariza- 
tion scheme prescribed by R.l, R.2 and R.3, we have shown the following: 
For both Case 1 and Case 2, the energy-momentum density T (x) has the 
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delta function expression (4-10), which reduces to T M (x) = -Mc 2 5°5 v 5^\x), 
when |0| = 1. 
Case 1 

— V 

(2.1. a) The correct expression for T ^ (x) can be obtained by transforming first 
the regularized density T „ (x; e) and then taking the limit e — > 0. However, 
it cannot be obtained if the order of the coordinate transformation and the 
limiting procedure is exchanged. 

(2.1.b) The regularized scalars R({}; e), RT({Y, e)R p „{{}; e) and R pa ^({}; e) 
xRpanu({}', e) are essentially equal to the corresponding scalars in Schwarzschild 
coordinates. 

Case 2 

(2. 2. a) The correct expression for T ^ (x) cannot be obtained by considering 
the transformation of the regularized or limiting densities in Schwarzschild 
coordinates. 

(2.2.b) The regularized scalars R({}; e), R» u {{}; e)i?^({}; e) and R papu ({}; e) 
xRpani>({}', e) are given by Eq. (4-24), from which Eqs. (4-25) and (4-26) 
follow. 

The following is worth mentioning: 

[A] When the condition \Q\ = 1 is satisfied, the gravity appears to be interpreted as 
produced by a point-like particle of mass M located at x = for both Case 1 and 
Case 2, so far as energy- momentum densities are concerned. Also for the orbit of a 
test point-like particle moving in the region r > a, this interpretation works well for 
Case 1 if the additional condition^] 

m(i-3<i> (5-i) 

following from Eq. (4-3), is satisfied. But, for Case 2, a test point-like particle is 
expected to move as if the source were located at r = —a/A. The fact that energy- 
momentum density T "(x) takes the form (4-11) does not necessarily imply that the 
gravity is produced by a mass point located at x = 0. 
10 Note that the condition ( |5-l| ) is satisfied for Schwarzschild and Kerr-Schild coordinates. 
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The situation is simple and natural in Case 1, as in Schwarzschild coordinates. This 
is reasonable, because the Jacobian of the transformation ( |1-4| ) for this case is given 
by 

QA 3 , (5-2) 



d(x) 

which implies that the transformation is everywhere regular. 

The volume element y /^gdx°dx 1 dx 2 dx 3 is given by 

^^dx°dx l dx 2 dx 3 = \ttA 3 \dx°dx l dx 2 dx 3 , (5-3) 

for Case 1, which includes Schwarzschild coordinates and Kerr-Schild coordinates as 
special cases. This case describes volume-preserving coordinates, which are preferred 
in the sense mentioned in Ref. [[J. 

For Case 2, the situation is rather curious, as is known from the statement [A]. 
The following is worth noting in this connection: The Jacobian of the coordinate 
transformation ( |1-4|) for this case is 

0{x) r z 

which means that the coordinate transformation is singular at r = and at Ar 
+a = 0. Thus, the coordinate systems of Case 2 and of Schwarzschild are not 
equivalent, and the corresponding metrics describe different physical situations. 6 -* 

The singularity at x = in Case 2 is weaker than in the case of the Schwarzschild 
coordinate system and in Case 1, as is seen from Eqs. (3-10), (3-11), (3-13), (4-18), 
(4-25) and (4-26). This is expected from the regularity and singularity of the Jaco- 
bians (5-2) and (5-4). 

We must be careful in considering the coordinate transformation of distributional 
quantities, as is known from Eqs. (3-9), (3-10), (3-11), (4-24), (4-25) and the state- 
ments (2.1. a) and (2. 2. a). 

The regularization schemes employed in §§ 3 and 4 are crucial in giving our results. We 
have regularized the "dynamical part" only. If we had regularized the "kinematical 



factors" 1/f and 1/f 2 in Eq. (3-4) as \j\/r 2 + e 2 and l/(f 2 + e 2 ), respectively, in 
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addition to regularizing the "dynamical function" h = a/f, for example, we would 
not have obtained Eq. ( |1-3| ). Also, if we had not followed the scheme prescribed by 
R.l, R.2 and R.3 in §4, we might have been led to unwelcome results. However, as 
for this scheme, we cannot simply say that we have regularized "dynamical parts" 
only. At present, we cannot give a lucid interpretation to the scheme employed in §4. 

Sometimes, a singularity is considered as a point existing outside of space-time and 
is thus removed from consideration. However, this way of handling such entities is 
not appropriate, because then the Schwarzschild gravity described by Eq. (1-1), for 
example, has to be regarded as produced by zero energy- momentum. A singularity 
exists in the space-time manifold, although it is not a "place". We believe it con- 
structive to examine space-time singularities by extending the analysis on manifolds 
in a distribution theoretical way, although we are not equipped with satisfactory 
mathematical machinery as yet. The expressions for scalar concomitants of the cur- 
vature tensor imply that mathematical formalism capable of incorporating objects 
like [S^ (x)~\ 2 and 5( 3 )(a;)/r 3 is needed for a precise description of singularities. 

We could not construct well-defined expressions for the quadratic scalars R pv {{}) 
xi? M „({}) and R paflu ({})Rpo-tiv({})- Otherwise, however, our prescription leads to 
reasonable results, such as Eqs. (1-3) and (3-10). Thus, the claim, "the well-known 
expression R patliy ({}) R pa ^ u ({}) = 12a 2 /f 6 is not correct • ■ ■" in the statement (l.B), 
is well-grounded. 

Both Case 1 and Case 2 do not satisfy the harmonic condition, and Case 1 with 
A = 1 = \Q\ is suited for the description of the Schwarzschild gravity and its source. 
This is opposed to the argument in favor of the harmonic condition as a physical 
supplement to Einstein's theory of gravitation. || |S], |Tt], [II]] For Lanczos coordinates, 
which satisfy the harmonic condition, the energy-momentum density of the source 
does not take the simple form Eq. (4-10). Harmonic coordinates are suited to discuss 
the gravitational wave, but they are not suited for the description of the Schwarzschild 
space-time. 
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Appendix A 

Proofs of relations (3-8) and (4-20) 

First, we prove the relation (3-8). Let us denneQ 



F a/3 (x; e) 



def X a X f3 



3 e 2 



+ 



2 [r 2 + e 2 ) 5 / 2 r 2 ( r 2 _|_ e 2 )3/2 
and let $(a?) be an arbitrary function of class C°° with compact support: 



p2w 



r 2 $ al3 (r) = I sin GdO I x a x p ®{x)d(j) = , for r > R 



to Jo 
with R being a positive constant. Then, we have 



J°0( e ) ^ y F Q/3 (a;;e)$(a;)d 3 a: = J ' *" '(r.r! 



2 (x 2 + l) 3 / 2 2 (a; 2 + If I 2 



which is calculated to give 



ra-1 



/^(e) = £ 



z=o 



Z! 



4 ^fc V 2 



1 + / Z 
2 



5 B „ 2 fi±i, 2 -i 
4 wr? V 2 2 



,n /■ — 



n\ ./o 



1 



_2(a; 2 + l) 3 / 2 2(x 2 + l) 5 / 2 _ 



dx , 



where we have expressed the function Q al3 (ex) as 

$^(ex) = V „ W (ecV + -±-(ex)"$ a/3(n) (0 , 
l\ n\ 

1=0 

i = 9ex , 1 > 9 > , n > 1 



(A-l) 



(A-2) 



dx , (A-3) 



(A-4) 



(A-5) 



with cjy^O) = f d l Q al3 /dr l . Also, -B 2 (p, g) represents the incomplete beta function of the first 
kind: 

(A-6) 
(A-7) 



B z (p,q) = / t p ~\l - t) q - l dt , 1 > Rez > 



Equation (A-4) gives 



\iml at3 (e) = 2tt5 q ^(x = 0) , 



11 In this Appendix, the symbol is omitted for simplicity. 
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where use is made of the relation 

$°0(O) = \n5 a ^{x = 0) . (A-8) 

Thus, we see that lim e ^ F a ^(x; e) = 27r5 a/ V 3 )(a;). This is equivalent to Eq. (3-8). 
Next, we prove Eq. (4- 20). Let us define 

G ^ e ) = VpT^' (A ' 9) 

and let ty(x) be an arbitrary function of class C°° with compact support: 

r m a (r) = / sin 6d0 / x a m(x)d<p = , for r > R > . (A-10) 
Jo Jo 

Then, we have 

//•— 2 
G^a; c)*(a ! )rf , a! = jf E ^ a (ex) | ds , (A-ll) 

from which the expression 

+ / <E«(")(fl dx, 9R>Z>0, (A-12) 

n\ Jo (x 2 + If I 2 

follows, where we have expressed the function ^> a (ex) as 

n_1 \b a ( l )(f)\ 1 

* a (^) = y^-^-^(exy + ^(exr^ n \o, 

A — ' 11 n\ 

1=0 

i = 9ex , 1 > 6 > , n > 1 (A-13) 

with = d l ^ a /dr l . Equation (A-12) leads to \im e ^ J a (e) = ^ a (0)/3 = 0, which gives 
Eq. (4-20). 

We note that the definite relation 



lim ± ^T FaP ( ay i e ) = 6vr5 (3) (a0 , (A-14) 

exists, while, the limit 



a, (3=1 



E -^H^^^ (A-15) 

a,/3=l 
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is not well-defined, as is known from the footnote |6] on page |9|. Also, we have 

3 r\ a 3 



a=l 



e^o z — ' r o z — ' r e— >0 



a=l 



X 

r 



We must be careful in treating F a/3 (x; e) and G a (x; e). 

Appendix B 

Expressions for the quantities t " , T ^ ' , R{{}), i? MZ/ ({})-R M „({}) and 



RW ({})R p<JIXV ({}) in terms of A,B,C and D 



(A-16) 



The energy-momentum density t is expressed as 



2/ct„ 



5 M , 2nt c 



BD 



2B'D' (B' 
BD + \ B 



1 A' 

rA 



r 



AB f 1 AC + £> 2 A' 

) ,J H < ' 

a VA{2r AB A 



A'B' 
! AB 



OC^ 1 '30 



where G is given by 



G 



AB 



A'B' 
~AB 



1 (B' 

2 V IF 



5' 
~B 



1 AC + D 2 A' 
r AB A" 



and A = A(B + C) + D 2 . Also, for T " , we have 



AB 



1 / A' B' A' 
r 1 A~ + 3 iF~ A 



AB 



+ 



1 / B' 



1 A'B' 1 B'A' 



kT 



0, «T„ 



AB 



4 VB 
o x Q BB> 



2 AB 2 BA 



r ^/A 



/ B' 




B" 




i 2 B 


40 


+ ~B 





IB' A' 
2 BA 



1 /A' B' 1A' 
r V A" + B" ~ 2 A" 



1A" IB" 1 /B' 
27 + 2 IF ~ 4 1 B" 



fB-r 



(B-2) 
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+ 



x a x f3 AB 



lA'B' 



2 AB 
2r A~ 



A' B'\ A' 
~A + ~b) A" 



1 AC + D 2 



1 /£' 
+ 2 B 



AB 
1 



A' 
7 



ia^ _ izr 

B'\ A'" 
~ B~J "a 



(B-3) 



For R({}),R^({})R^({}) and i?^({})i?^({}), we have 



R({}) = £ 



2 

r 

+: 



A' B' A' 



+ 



1 (B[_ 
~B 



MB' 
2 AB + 



2 AC + D 2 _ 
AB 

B'\ A' 
A 



A" B" 
~A~ ~B 



1M 

2A + ~B 



A 2 /I A" 



1 A' A' 1 A'B' 1 A' 
4 AA + 2^8" + rl 



A 2 
+2 A~ 2 



A 2 
+ A^ 



1/1 A' A' B' 
r \2~K ~A~ 2 ~B 



1 B" 1 B'A' 
W + 4 BA 



+ 



1 AC + B 2 1 A'B' 



AB 



2 AB 



I A" IMA' 1MB' B" IfB' 
2T~4^4A" + 2i5 + 2 I IF 



i^({})W{}) = 



A 2 (M_ 
A 2 

A 2 
A2 

A 2 
A 2 



1 5' A' 1 
"2 BA + r 



1 A'A' 



A' A' 5' 
^A^B 



A 2 



+4 
+2 



2 AA 
15' 1 AC + D 2 



. ,1M 1MB' 

+ 2 — 1 

A 2 V r A 2 AB 



r B 

1 

r 

1 

"2 



r 2 AB 

A' B' A' 
1 + 2 l3- A 



+ 



B' 
IF 



1 A'B' B" 
+ 2 + B" 



IF 



IB' A' 
2 BA 



(B-4) 



In the above, we have defined M = dA/dr, A" = d 2 A/dr 2 , etc. 
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